Spacetime foam as a quantum thermal bath 
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An effective model for the spacetime foam is constructed in terms of nonlocal interactions in a 
classical background. In the weak-coupling approximation, the evolution of the low-energy density 
matrix is determined by a master equation that predicts loss of quantum coherence. Moreover, 
spacetime foam can be described by a quantum thermal field that, apart from inducing loss of 
coherence, gives rise to effects such as gravitational Lamb and Stark shifts as well as quantum 
damping in the evolution of the low-energy observables. 
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The foamlike structure of spacetime was first suggested 
by Wheeler ]jj and, since then, various components, such 
as wormholes an d virtual black holes have been 
proposed. The quantum fluctuations of the geometry 
that constitute the spacetime foam should be of the same 
order as the geometry itself at the Planck scale. This 
would give rise to a minimum length Q beyond which 
the geometrical properties of spacetime would be lost, 
while on larger scales it would look smooth and with a 
well-defined metric structure. 

Planck length Z* might play a role analogous to the 
speed of light in special relativity. In this theory, there 
is no physics beyond this speed limit and its existence 
may be inferred through the relativistic corrections to the 
Newtonian behavior. This would mean that a quantum 
theory of gravity could be constructed only on "this side 
of Planck's border" as pointed out by Markov ||]. In fact, 
the analogy between quantum gravity and special relativ- 
ity is quite close: in the latter you can accelerate forever 
even though you will never reach the speed of light; in 
the former, given a coordinate frame, you can reduce the 
coordinate distance between two events as much as you 
want even though the proper distance between them will 
never decrease beyond Planck length (see Ref. fjj, and 
references therein). This uncertainty relation Ax > Z* 
also bears a close resemblance to the role of Ti in quan- 
tum mechanics: no matter which variables are used, it is 
not possible to have an action / smaller than h. Indeed, 
the uncertainty principle can adopt the form AI > h. 

Spacetime foam and the related lower bound to space- 
time uncertainties would leave their imprint in low- 
energy physics. Indeed, low-energy experiments would 
effectively suffer a nonvanishing uncertainty coming from 
this lack of resolution in spacetime measurements. Then 
a loss of quantum coherence would be almost unavoidable 
Qj. It could also be expected that other effects such as 
transition-frequency shifts and quantum damping, char- 



acteristic of systems in a quantum environment H , may 
be present. In this Letter, we in fact show that spacetime 
foam behaves as a quantum thermal bath with a nearly 
Planckian temperature. 

In order to build an effective theory, we will substitute 
the spacetime foam, in which we possibly have a mini- 
mum length because the notion of distance is not valid at 
such scale, by a fixed background with low-energy fields 
living on it. We will perform a 3+1 foliation of the ef- 
fective spacetime that, for simplicity, will be regarded 
as flat, t denoting the time parameter and x the spa- 
tial coordinates. The gravitational fluctuations and the 
minimum length present in the original spacetime foam 
will be modeled by means of nonlocal interactions that 
relate spacetime points that are sufficiently close in the 
effective background, where a well-defined notion of dis- 
tance exists. Furthermore, these nonlocal interactions 
will be described in terms of local interactions as follows. 
Let {/ii[t]} be a basis of local gauge-invariant interac- 
tions at the spacetime point [x, t) made out of factors 



of the form ll n ^ 1+s ^ 4 [<f>(x, t)Y' L , <j> being the low-energy 
field strength of spin s. As a notational convention, each 
index i implies a dependence on the spatial position x; 
also any contraction of indices will also entail an integral 
over spatial positions. Then, we can write the nonlo- 
cal effective interaction term in the Euclidean action as 
lint = J2n i n with 



N 



In = W\' dt i--- dt Nc il " 4N (t 1 ...t N )h il [ti]---hi N [t 
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Here, c n "' lN (t\ . . . t^) are dimensionless functions that 
vanish for relative spacetime distances larger than the 
length scale r of the gravitational fluctuations. Fur- 
thermore, these coefficients can depend only on rela- 
tive positions and not on the location of the gravita- 
tional fluctuation itself. The physical reason for this 
is conservation of energy and momentum: the fluctua- 
tions do not carry energy, momentum, or gauge charges. 
Thus, diffeomorphism invariance is preserved, at least at 
low-energy scales. One should not expect that at the 
Planck scale this invariance still holds. However, this vi- 
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olation of energy-momentum conservation is safely kept 
within Planck scale limits |l0|, where the processes will 
no longer be Markovian. Furthermore, the coefficients 
c n - iN (h . . .t N ) will also contain a factor [e~ s ^/ 2 ] N , 
S(r) being the Euclidean action of the gravitational fluc- 
tuation, which is of the order (r/l*) 2 . 

Since higher-spin s > or higher-power n > 1 interac- 
tions are suppressed by inverse powers of the low-energy 
length scale I, we will concentrate on the mass term for 
scalar fields hi[t] — l^^x 1 , t) 2 , where now the index i 
just keeps track of the dependence on the spatial posi- 
tion. A simple calculation shows that In ~ £ N (l/ r ) 4 , 
where e = e -5 '( , ')/ 2 (r/Z >t ) 4 (Z/Z»)~ 2 . The parameter e has 
contributions from two different scales: on the one hand, 
it depends on the length scale of the gravitational fluctu- 
ations r and, because of the exponential factor, it will be 
very small for fluctuations of few Planck lengths; on the 
other hand, it depends on the low-energy scale I through 
its inverse squared and will therefore be very small far 
from Planck's regime. In the weak-coupling approxima- 
tion, i.e., up to second order in the expansion parameter 
e, the trilocal and higher effective interactions do not 
contribute. The terms To and I\ are local and can be 
absorbed in the bare action (note that the cofficient c 
appearing in 2"o is constant and that the coefficients d (t) 
in X\ cannot depend on spacetime positions because of 
diffeomorphism invariance). Consequently, we can write 
the interaction term as a bilocal term in the Euclidean 
action 

hr* = \ J dtdt' c«(t-t')MWL 

where c y (f — t') is of order e~ s ^ and is concentrated 
within a spacetime region of size r. Then, the effective 
partition function has the form Z — J T>4> e _/o+7int , I Q 
being the bare low-energy action for the scalar field. 

This bilocal effective action, when rotated back to 
Lorentzian spacetime, does not lead to a unitary evo- 
lution. The reason for this is that it is not sufficient to 
know the fields and their time derivatives at an instant 
of time in order to know their values at a later time: 
we need to know the history of the system. There exist 
different trajectories that arrive at a given configuration 
(</>, <fi). The future evolution depends on these past tra- 
jectories and not only on the values of </> and (j) at that 
instant of time. Therefore, the system cannot possess a 
well-defined Hamiltonian vector field and suffers from an 
intrinsic loss of predictability |Tf| . This can be best dealt 
with by writing, up to a determinant, the exponential of 
the interaction term as jl2| 

glint ^ J D ae -\ I dtdt'njit-t')^ (t) a j (O e - J dta* (t)hi[t]_ 

Here, the continuous matrix "fij(t — t') is the inverse of 
c«(t - t'), i.e., / dt" llk {t - t")c k i{t" - t') = 5(6(1 - t 1 ). 



We see that a is a random spacetime function subject to 
a Gaussian distribution. At second order in e and lowest 
order in r/l, the two-point correlation function is equal to 
{a l {t)ai(t')) = c ij (t - t') and (a 1 it)) = 0. Note that the 
Gaussian character of the distribution for the noise a is 
a consequence of the weak-coupling approximation (sec- 
ond order in e) , which keeps only the bilocal term in the 
action. Higher-order terms would introduce deviations 
from this noise distribution. The nonunitary nature of 
the bilocal interaction has been encoded inside the func- 
tion a, so that, when insisting on writing the system in 
terms of a Hamiltonian, an additional sum over the part 
of the system that is unknown naturally appears. Note 
also that we have a single function a 1 it) because we are 
considering only one local interaction; we will have a dif- 
ferent function a for each kind of interaction. 

The Lorentzian dynamics of the low-energy field will be 
governed by a master equation which can be derived af- 
ter a number of steps and approximations that are briefly 
outlined in what follows. For each fixed function a, we 
first calculate the evolution equation for the density ma- 
trix p a {t) obtained with the Hamiltonian 

H a (t) = H Q [t]+a l (t)h t [t], 

Hq [t] being the bare Hamiltonian of the low-energy field, 
and transform this equation into the interaction picture. 
We then integrate this equation between and t with two 
iterations and differentiate the result, so that the evolu- 
tion equation becomes an integro-differential equation for 
the density matrix. Next, we perform the Gaussian av- 
erage over a and expand the result up to second order 
in the parameter e, taking into account that p a does not 
depend on a at zeroth order but only at first order in 
e, i.e., p a — p + 0(e) with p = (p a ) (weak-coupling ap- 
proximation). We also assume that pit) hardly changes 
within a correlation time r (Markov approximation), so 
that pit + r) ~ pit)- Finally, we transform the resulting 
equation back to the Schrodinger picture. At the low- 
est order in r/l, the result is a master equation for the 
low-energy density matrix which has the form Jf3| 

/>oo 

p = -i[H ,p] - / dTC ij (r) [hi, [hj,p\] . 
Jo 

The first term gives the Hamiltonian evolution that 
would also be present in the absence of fluctuations. The 
second term is a diffusion term which will be responsible 
for the loss of coherence (and the subsequent increase of 
entropy) . It is a direct consequence of the foamlike struc- 
ture of spacetime and the related existence of a minimum 
length. 

The characteristic decoherence time induced by the 
diffusion term can be easily calculated and yields the fol- 
lowing ratio between the decoherence time and the low- 
energy length scale: t&/1 ~ e s ^ (r /h.)~ 4 ■ Because of 
the exponential factor, only the gravitational fluctuations 
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whose size is very close to Planck length will give a suffi- 
ciently small coherence time. Slightly larger fluctuations 
will have a very small effect on the unitarity of the ef- 
fective theory. For higher spins and/or powers of the 
field strength, the decoherence time increases by powers 
of I /l*. For instance, if we consider interactions that mix 
fields with different spin, then the next relevant decoher- 
ence time corresponds to the scalar-fermion interaction 
term 4> 2 4>ip, which has an associated decoherence ratio 
Td/l proportional to I /I*. Note that this decoherence time 
may be small enough for sufficiently high energies. 

Let us now go a bit further and describe spacetime 
foam in terms of a quantum thermal bath. With this aim, 
we will consider a system consisting of the low-energy 
fields coupled to a quantum bath || . By comparing this 
system with the results obtained above for gravitational 
fluctuations, we will see that the latter can be substituted 
by a thermal bath. So, let us start with a Hamiltonian 
of the form 

H = Hq + H int + H[,. 

Hq is the bare Hamiltonian that represents the low- 
energy fields and Hb is the Hamiltonian of a bath that, 
for simplicity, will be represented by a massless scalar 
field. The interaction Hamiltonian will be of the form 
■Hint = C^ii where the noise operator £ is of the form 
f (t) = i J dk^ X {u)[a + {k)e l ^ l - kx ^ - H.c.]. In this ex- 
pression, a and a + are, respectively, the annihilation and 
creation operators associated with the bath, u> — V«?, 
and x(w) represents the coupling between the system 
and the bath. This implies that £'(t) = with 
Pj (t) = p(x° , t) being the momentum of the bath scalar 
field and x u = J dkx{uj) cos[k(x l — being the cou- 
pling between the low-energy field and the bath in the 
position representation. The coupling x( w ) must be such 
that there exists a significant interaction with all the bath 
frequencies u up to the natural cuttoff r -1 . All the rel- 
evant information about the coupling is encoded in the 
commutation relations and the correlation function of the 
noise operator £. 

Since the commutator of the noise operator £ at dif- 
ferent times is a c number, we can introduce the so- 
called commutative noise representation @, which will 
allow us to compare this model with that of topological 
fluctuations previously described. This can be done by 
defining a new noise operator a in the following form: 
a l {t)p{t') ee i[f (t),p(t')]+. It is straightforward to 
check that the operator a commutes at any time, i.e. 

a? (t')] — 0. However, this does not mean that it 
commutes with everything. Indeed, the commutator of a 
with any low-energy operator A is nonvanishing and has 
the form: 

[A\t),ai(t')} = [ dT[A i {t),h k (r)]f ik V-T), 
Jo 

where 



f ij (r) = / du>u 2 G ij (w) cos ujt, 
Jo 



The commutator above vanishes for low-energy opera- 
tors that are in the far past of the noise and is nonzero 
when they are in the near past or the future. Only in the 
so-called first Markov approximation the frontier among 
both regimes is sharply located where both noise and 
low-energy fields are at the same instant of time. There- 
fore, the function /^(t) can be interpreted as a kind of 
memory function. 

If we assume that the bath is in a thermal state pb = 
Z~ 1 e~ Hb ' T with a temperature T and define the average 
of any operator Q as (Q) = Tvb(Qpb), we can compute 
the correlation function c* J '(f — t') ee (a % (t)a? (t'))\ 

poo 

c 13 (t) = / duju 3 G ll {uj)[N(uj) + 1/2] coswr, 
Jo 

where N(uj) = [exp(w/T) — 1] 1 is the mean occupation 
number of the bath corresponding to the frequency u>. 
Also, it can be shown that the trace (Q) corresponds to 
a Gaussian average over a only in the case that the bath 
is in a thermal state || , as we are considering. 

We are now ready, following similar steps to those out- 
lined before, to write down the master equation for the 
low-energy density matrix. If we keep terms only up to 
second order in the expansion parameter e given by the 
product of the thermal correlation time r a ~ 1/T of a, 
the size of the operator h and the root mean square of 
a, which is of the order of y/B, and we also assume that 
T & l I -C 1, then the resulting equation has the same form 
as the classical master equation obtained above with the 
correlation function c y (r) substituted for c y (r). From 
the identification of both models (a = a), we conclude 
that the temperature of the heat bath is determined by 
the size of the gravitational fluctuations, i.e., T ~ 1/r 
and that e = e <C 1 (weak coupling approximation). 
Note also that the coupling x( w ) is uniquely determined 
by the correlation function c y (r) by means of a suitable 
mode expansion. The zeroth order approximation hxr/l 
that we have made in order to compare and identify both 
models can be regarded as a kind of semiclassical approx- 
imation since all the quantum features of the noise have 
disappeared from the master equation. 

We can however obtain a more general master equa- 
tion, valid up to second order in e and with no restriction 
bxr/l, that takes into account the quantum nature of the 
gravitational fluctuations. These contributions will be 
fairly small in the low-energy regime, but may provide 
interesting information about the higher-energy regimes 
in which I may be of the order of a few Planck lengths 
and for which the weak coupling approximation is still 
valid. The quantum noise effects || are reflected in the 
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master equation through a term proportional to f lJ (r) 
and another proportional to (?i(r), both of them inte- 
grated over r 6 [0, oo]. Because of these incomplete inte- 
grals, each term provides two different kinds of contribu- 
tions whose origin can be traced back to the well-know 
formula J Q d,Te lu}T = it8(lu) + P(i/uj), where P is the 
Cauchy principal part [jlij. Thus, the / term contains a 
dissipation part, necessary for the preservation of com- 
mutators, and a contribution to what can be interpreted 
as a gravitational Lamb shift. On the other hand, the 
c-term gives rise to four different contributions: The al- 
ready discussed diffusion term, another diffusion term 
originated from the vacuum fluctuations of the bath and 
that does not vanish at zero temperature, another con- 
tribution to the gravitational Lamb shift and, finally, a 
shift in the scalar-field oscillation frequencies that can be 
interpreted as a gravitational Stark effect. The size of 
these effects, compared with the bare evolution, can be 
calculated after some work: the thermal diffusion term is 
of order e~ s ^ (r /h ) 4 , which is the only one that survived 
in the previous approximations; the diffusion created by 
vacuum fluctuations, the damping term, and the Stark 
effect are smaller by a factor r/l; and the Lamb shift is 
smaller than the diffusion term by a factor (r/l) 2 . 

The models described in this Letter are particularly 
suited to the study of low-energy effects produced by 
simply connected topology fluctuations such as virtual 
black holes Q. Indeed, it has been shown that virtual 
black holes can be represented from the low-energy point 
of view by effective interactions h{[t] like the ones em- 
ployed here. The master equation can then be interpreted 
as providing the evolution of the density matrix in the 
presence of a bath of ubiquituous quantum topological 
fluctuations of the virtual-black-hole type. 

Multiply connected fluctuations (with vanishing sec- 
ond Betti number) such as wormholes (|] can also be 
described as nonlocal interactions that, in the weak- 
coupling approximation, become bilocal. The coefficients 
c lJ of this bilocal term do not depend on spacetime posi- 
tions since multiply connected topology fluctuations con- 
nect spacetime points that may be far apart from each 
other. Diffeomorphism invariance also requires the space- 
time independence of c lJ . This can also be seen by an- 
alyzing these wormholes from the point of view of the 
universal covering manifold, which is by definition sim- 
ply connected. Here, each wormhole is represented by 
two boundaries located at infinity and suitably identi- 
fied. This identification is equivalent to introducing co- 
efficients c ij that relate the bases of the Hilbert space 
of wormholes in both regions of the universal covering 
manifold. Since c 4 -? are just the coefficients in a change 
of basis, they will be constant. As a direct consequence, 
the correlation time for the functions a 1 is infinite. This 
means that the functions a 1 cannot be interpreted as 
noise sources that are Gaussian distributed at each space- 
time point independently. Rather, they are infinitely co- 



herent. The Gaussian distribution to which they are sub- 
ject is therefore global, spacetime independent [||. Con- 
sequently, the master equation contains no diffusion term 
and, actually, it predicts a unitary evolution for the den- 
sity matrix. If we still try to represent wormholes by 
a thermal bath as we have done with localized gravi- 
tational fluctuations, we soon realize that, in order to 
reproduce the infinite correlation time, the couplings 
must be constant, that they must commute with every 
other operator and, related to these two facts, that only 
the zero-frequency mode of the bath can be coupled to 
the low-energy fields, in agreement with the result that 
the Gaussian distribution is spacetime independent and 
that the effective theory is, in this case, unitary. 

Let us conclude with a brief summary. We have de- 
scribed spacetime foam in terms of a quantum thermal 
field, which induces a loss of coherence in the low-energy 
dynamics as well as other effects of quantum nature such 
as dissipation and gravitational Lamb and Stark shifts. 
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